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ABSTRACT 

The orbits of the outer two known planets orbiting Upsilon Andromedae are 
remarkably eccentric. Planet C possesses an orbital eccentricity of e\ = 0.253. 
For the more distant planet D, €2 = 0.308. Previous dynamical analyses strongly 
suggest that the two orbits are nearly co-planar and are trapped in an apsidal 
resonance in which Auj, the difference between their longitudes of periastron, 
undergoes a bounded oscillation about 0°. Here we elucidate the origin of these 
large eccentricities and of the apsidal alignment. Resonant interactions between 
a remnant circumstellar disk of gas lying exterior to the orbits of both planets 
can smoothly grow ei- Secular interactions between planets D and C can siphon 
off the eccentricity of the former to grow that of the latter. Externally amplifying 
€2 during the phase of the apsidal oscillation when e^/e-i is smallest drives the 
apsidal oscillation amplitude towards zero. Thus, the substantial eccentricity of 
planet C and the locking of orbital apsides are both consequences of externally 
pumping the eccentricity of planet D over timescales exceeding apsidal precession 
periods of order 10 4 yr. We explain why the recently detected stellar companion 
to v And is largely dynamically decoupled from the planetary system. 

Subject headings: planetary systems — celestial mechanics — stars: individual 
(v Andromedae, HD168443, 47 UMa) 



1. INTRODUCTION 

The outer two known planets orbiting Upsilon Andromedae (v And) possess remarkably 
large orbital eccentricities (Butler et al. 1999). Planet C has a mass of mi = (1.83/ sinii) Mj 
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(where %\ is the angle between the orbit pole of planet C and our line of sight, and Mj is the 
mass of Jupiter), an orbital semi-major axis of a\ = 0.805 AU, and an orbital eccentricity 
of ei = 0.253. For planet D, the corresponding quantities are m 2 = (3.79/ sin i 2 ) Mj, a 2 = 
2.48 AU, and e 2 = 0.308. These orbital parameters were kindly supplied by D. Fischer 
(2002, personal communication), and represent more up-to-date values than those employed 
by previous works. Studies of the dynamical stability of the system (Rivera & Lissauer 2000; 
Stepinski, Malhotra, & Black 2000; Lissauer & Rivera 2001; Chiang, Tabachnik, & Tremaine 
2001) strongly suggest that the two planets execute approximately co-planar trajectories that 
are observed nearly edge-on, so that sinii pa shu 2 > 0.5. 

That m 2 likely exceeds m x while e 2 > e 1 stands at odds with the idea that gravitational 
interactions between planets D and C excited both e 2 and e x to their present-day values. 
Close encounters between planets (Rasio & Ford 1996; Weidenschilling & Marzari 1996; 
Ford, Havlickova, & Rasio 2001) and mean- motion resonance crossings between planets on 
divergent orbits (Chiang, Fischer, & Thommes 2002) tend to impart the greater eccentricity 
to the less massive body. Thus, we are led to the conclusion that an external agent— 
another planet, a star, or the circumstellar disk from which the planets formed — must have 
played a direct role in exciting the eccentricity of the most massive planet, D. Of these three 
candidates, the latter is the most promising, as we explain below. 

An additional clue to the origin of v And lies in the apsidal resonance that obtains if 
«2i, the mutual inclination between the orbits of planets C and D, is less than ~20°. 3 The 
observation today that the arguments of periastron of the two orbits, ou\ and u 2 , differ by 
only ~12° has been used by Chiang, Tabachnik, & Tremaine (2001, hereafter CTT) to argue 
that the orbits are indeed nearly co-planar. 4 If i 21 < 20°, then Alu = uo 2 — uj 1 librates about 
0° with a semi-amplitude of ~38°. This apsidal libration and the eccentricity evolution of 
both planets are displayed in Figure 1. 

This paper lays the groundwork for understanding the origin of the large eccentricities 
and the apsidal alignment exhibited by the orbits of planets C and D in v And. Our 
main result is that the eccentricity of planet C and the locking of orbital apsides are both 
consequences of the slow growth of the eccentricity of planet D. The latter eccentricity, 
in turn, was driven by an external agent — plausibly a primordial circumstellar disk lying 



3 The mutual inclination, i 2 i, does not necessarily equal i 2 — h; see, e.g., Stepinski, Malhotra, & Black 
(2000). 

4 Thc argument of periastron, u>, as fitted by Doppler velocity measurements is referred to the plane of 
the sky (see, e.g., CTT). The longitude of periastron, Q, may be referred to any fixed line in inertial space 
for the purposes of this paper. When i 2 i =0°, we can take ui = u. 
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Fig. 1. — Time evolution of orbital eccentricities of planets C (denoted 1) and D (denoted 
2), and the difference in their apsidal longitudes, as computed using linear secular theory 
for sinix = sini 2 = 1 and i 2 i — 0°. The maximum in ou 2 — &i is 38°; this differs from 
the preferred value cited in CTT because of the inclusion of more recent Doppler velocity 
measurements. The system spends more time at small e2/ei than at large e-ije^ a fact that 
is important for understanding the damping of apsidal libration. 
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exterior to the orbit of planet D — that acted over timescales exceeding 10 4 yr. We play 
our scenario out and explain the mechanics of apsidal resonance capture in §2. We also 
discuss in that section why the recently detected companion star to v And is not likely to 
be dynamically relevant. A summary is provided in §3. 

2. Secular Growth of e\ and Apsidal Resonance Capture 

2.1. Basic Model 

We take planets C and D to occupy co-planar, nearly circular orbits at time t = 0. The 
dynamics exhibited by the system is qualitatively the same over the range 0° < i 2 i < 30° 
(CTT). We adopt current observed values for the orbital semi-major axes, and minimum 
values for the planetary masses (sin^ = sini 2 = 1)- 

We introduce an external force only on planet D that smoothly amplifies that planet's 
eccentricity over a finite time interval. The equations governing the eccentricities and apsidal 
longitudes (a)) of the orbits of planets C and D read as follows (cf. Murray & Dermott 1999): 



hi EE 


|(e lS in^) 


= +A 11 k 1 + A l2 k 2 


(1) 


k\ EE 


— (ei cosc^i) 


= -Auh! - A 12 h 2 


(2) 


h 2 = 


|(e 2 sin^) 


= +A 21 k t + A 22 k 2 + Eh 2 


(3) 


k 2 = 


|(e 2 cos^ 2 ) 


= —A 21 h\ - A 22 h 2 + Ek 2 


(4) 



where 



A 11 
A 12 
A 21 
A 22 



Hi 

4 + mi 
rii m 2 
4 m t + mi 

n 2 mi 
4 m* + m 2 
n 2 mi 



4 m* + m 2 



J 3/2 



6 (2) 
U 3/2 



~ 1 b 3/2 I — 



J 3/2 



(5) 
(6) 
(7) 
(8) 



- 5 - 



E = exp(-t/r E )/r e . (9) 

Here rij is the mean motion of planet j, m* = 1.3M is the central stellar mass, b^j 2 and bf^ 2 
are the usual Laplace coefficients, and 1/E is the timescale over which exponential growth 
of e<i occurs, where te and r e are fixed time constants. 

We emphasize that only the eccentricity of planet D is directly amplified by our external 
force. Our equations are intended to represent the following physical picture, staged just 
after the formation epoch of the planets. Interior to and between the orbits of planets C and 
D, circumstellar disk gas is absent. Only exterior to the orbit of planet D lies disk gas whose 
resonant interaction with that planet excites the planet's eccentricity. Our scenario is con- 
sonant with most current thinking on planet-disk interactions. Gas between the two planets 
is plausibly driven out of this region by planet-induced torques, as numerical simulations by 
Kley (2000) and Bryden et al. (2000) find. Such gas either tunnels past both planets into the 
innermost and outermost disks, or is accreted by the planets (Bryden et al. 2000). Viscous 
accretion may drain away the disk that lies interior to both planets onto the central star, 
as has been invoked by Snellgrove, Papaloizou, & Nelson (2001). The disk outside ~3 AU 
may be prevented from accreting inwards by shepherding torques exerted by the outermost 
planet. Goldreich & Sari (2002) have discovered that disk gas can excite a planet's eccentric- 
ity; eccentricity-amplifying Lindblad resonances can defeat eccentricity-damping co-rotation 
resonances provided the planet's initial eccentricity exceeds a critical threshold. While the 
nonlinear outcome of their finite amplitude instability remains to be worked out, it may be 
that a planet's eccentricity ceases to grow when the planet's orbit begins to overrun the gap 
edge, i.e., when e ~ 0.3. For a dimensionless disk viscosity of a = 10~ 4 and a disk mass of 
40Mj, we compute an eccentricity amplification timescale, e-ije^ of 7 x 10 4 yr, and a radial 
migration timescale, Ic^/a^l; of 4 x 10 6 yr [Goldreich & Sari 2002; their equations (18) and 
(19)]. Since the latter can exceed the former, we fix the semi-major axes of both planets in 
our analysis. 5 

We are aware that certain details underlying our scenario are not well understood. Gas 
between planets C and D might be shepherded rather than driven out (Bryden et al. 2000), 
in which case its effects on D and C would need to be accounted for. Numerical simulations 
by Papaloizou, Nelson, & Masset (2000) of interactions between Jupiter-mass planets and 
their parent disks result in eccentricity damping and have yet to be reconciled with the 
predictions of Goldreich & Sari (2002). Uncertainties in the plausibility of our scenario 
notwithstanding, the key features of equations (l)-(4) are that (a) the eccentricity evolution 



5 We have verified by direct numerical simulations that relaxing this assumption changes none of our 
conclusions. 
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of planet D is directly driven by an external influence, and (b) the eccentricity and apsidal 
evolution of planet C are not, i.e., planet C feels only the gravitational potential of planet 
D. These minimal requirements are sufficient to understand the origin of the eccentricity of 
planet C and the origin of the curious apsidal alignment of the orbits of C and D, as we now 
demonstrate. 

If E = 0, equations (l)-(4) yield the classical Laplace-Lagrange (L-L) solution. One 
solution, fitted with initial conditions appropriate to observations of v And today, is em- 
bodied in Figure 1. A variety of L-L solutions, differing only in initial eccentricities and 
apsidal longitudes, are mapped in the space of Ad; = Co 2 — Cj\ and e 2 /ei in Figure 2a. The 
dashed separatrix divides circulating solutions from librating solutions. Librating solutions 
all have \Auj\ < 90°. Today planets C and D live on a contour of libration characterized by 
max I Aa>| m 38°. 

If 0(E) <C O(Ajk), or equivalently if 1/E is longer than apsidal precession timescales of 
order 10 4 yr, then the system evolves adiabatically through a series of L-L solutions. Figure 
2b follows the evolution of a system for which ei(t = 0) = e 2 (t = 0) = 0.05, Au>(t = 
0) = 180°, r e = 7.0 x 10 4 yr, and t e = 1.5 x 10 5 yr. The final time of the calculation is 
tf — 1 x 10 6 yr, at which time the eccentricity driving term, E, is effectively zero. Each 
point represents the position of the system in Auj — e 2 j ' e\ space at a particular time; the 
solution was obtained by numerically integrating equations (l)-(4) using a standard fourth- 
order Runge-Kutta routine (Press et al. 1992). A set of points highlighted by a symbol of a 
given type indicates the trajectory over a time interval of 8400 yr. Comparison of Figure 2a 
with 2b indicates that the system morphs smoothly through a series of L-L contours. What 
is noteworthy is that the system is captured into libration and thereafter evolves towards 
smaller libration amplitude. 

The mechanics underlying apsidal resonance capture and the damping of apsidal libra- 
tion is as follows. The introduction of our eccentricity driving term amplifies e 2 , i.e., it 
pushes the system to the right in Figure 2. Thus, a system initially occupying a contour of 
circulation is eventually pushed past the separatrix onto a contour of libration. It remains 
to explain the evolution towards smaller libration amplitude. A librating system executes an 
elliptical trajectory in Ad; — e 2 /ei space. While the system traces out the left-hand side of 
an elliptical contour, external amplification of e 2 pushes the system to the right, i.e., towards 
contours of smaller libration amplitude. While the system traces out the right-hand side of 
the ellipse, external amplification of e 2 also pushes the system to the right, but now towards 
contours of larger libration amplitude. Which half of the interaction dominates depends on 
how much time the system spends on each half of the ellipse. It is clear from Figure 1 that 
the system spends more time near the minimum value of e 2 /ei than near the maximum 
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(a) Laplace-Lagrange Contours 
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(b) Computed Trajectory 




Fig. 2. — (a) Solutions of equations (l)-(4) for E = 0. Each contour represents one possible 
solution. The dashed contour closes at infinite e 2 /ei; it is the separatrix dividing apsidally 
circulating solutions from apsidally librating ones, (b) Solution of equations (l)-(4) for 
ei (t = 0) = e 2 (t = 0) = 0.05, ACu(t = 0) = 180°, r e = 7.0 x 10 4 yr, and r E = 1.5 x 10 5 yr, 
over a time span of 1 x 10 6 yr. Each point marks the instantaneous position of the system 
in phase space. Open circles denote the trajectory over a 8400 yr interval starting at t s = 0; 
+'s, t s = 5.1 x 10 4 yr; X's, t s = 1.2 x 10 5 yr; and filled circles, t s = 2.8 x 10 5 yr. The system 
is captured from circulation into libration and evolves towards smaller libration amplitude. 
Damping of libration results from increasing e 2 when e± is smallest [filled arrows in (a)]. 
Increasing e 2 when e 2 /ei is largest excites libration [open arrows in (a)], but this effect is 
small because the system spends less time at large e 2 /ei than at small e 2 /ei. 
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value of e 2 /ei. The same behavior can be seen from Figure 2b; the points, plotted at equal 
time intervals, cluster more strongly near small eije\ than near large e^jex. The reason for 
this is that the apsides of low eccentricity orbits tend to precess more quickly than those of 
high eccentricity orbits; for the same perturbative acceleration, oj oc — e 2 /e by Gauss's 
equation; and thus when e\ is minimal (e2/ei maximal), the orbital apsides of planet C 
precess more rapidly past the apsides of planet D than when e\ is maximal (e2/ei minimal). 
Variations in e 2 can be ignored in this analysis; they are characteristically smaller than vari- 
ations in ei because planet D carries the lion's share of the orbital angular momentum of 
the system (see, e.g., CTT and Figure 1). Thus, external driving of e 2 on the leftmost half 
of an elliptical contour dominates the driving of e 2 on the rightmost half, and the system 
spirals deeper into apsidal lock — all without the need for any explicit dissipation of energy! 

The celestial mechanical considerations described above are general and we have veri- 
fied by numerous integrations of equations (l)-(4) that apsidal resonance capture and the 
subsequent damping of apsidal libration occur under a wide variety of initial conditions. 
The input parameters employed for the trajectory shown in Figure 2b were chosen so that 
the end state of our calculation yielded system parameters that resemble those of v And 
today. Figure 3 portrays the same model, displaying the eccentricity and apsidal evolution 
as explicit functions of time, and should be compared with Figure 1. 

2.2. Refinements and Extensions 

2.2.1. Disk-induced Precession 

In addition to exciting the eccentricity of planet D, an external disk would also cause the 
orbits of both planets to precess at rates different from those due solely to the planets' mutual 
gravity. We may model disk-induced precession by adding extra terms to An and A 2 2 that 
decay on timescales characterizing the dissipation of the disk (by, e.g., photoevaporation). 
This procedure can account for the disk's temporary resonant and secular contributions to 
planetary precession rates. We have experimented by adding various transient terms to An 
and A 2 2 and find that none of our conclusions changes qualitatively. Faced with extra, disk- 
driven differential precession, we may still reproduce the observed parameters of v And by 
appropriately re-adjusting our input parameters (e.g., initial eccentricities and r e ). 
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Fig. 3. — Solution of equations (l)-(4) using input parameters listed under Figure 2. Exter- 
nally driving the eccentricity of planet D amplifies the eccentricity of planet C and locks the 
system into apsidal resonance. At the end of the integration, the eccentricities and apsidal 
longitude difference match those of v And today. 
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2.2.2. Stellar Companion 

A stellar companion to v And has recently been reported by 2MASS (Two Micron 
All-Sky Survey; Lowrance, Kirkpatrick, & Beichman 2002). The mass of the companion is 
m E ~ O.2M and the projected separation between it and the primary is 750 AU. Could 
this companion star be responsible for exciting the observed eccentricities of planets C and 
D? 

The answer is almost certainly no. The usual Kozai mechanism for pumping planetary 
eccentricities (see, e.g., Holman, Touma, & Tremaine 1997) can only operate if (1) the plane 
of the planetary system is sufficiently inclined with respect to the plane of the stellar binary, 
and (2) the apsidal precession rate of the planet due to the stellar companion dominates 
other contributions to that precession rate. The latter condition fails to be satisfied, by 
a wide margin. We define ce and oe to be the (unknown) orbital eccentricity and semi- 
major axis, respectively, of the companion star. Then the precession period of planet D, as 
induced solely by the stellar companion, is of order ~ 27m 2 ~ 1 (m* / m E ) (ag / a 2 ) 3 '(1 — e^) 3 / 2 ~ 
2.5 x 10 8 (0.2M o /m s )(a B /750 AU) 3 [(1 - e E )/0.5} 3/2 yr (Holman, Touma, & Tremaine 1997). 
This greatly exceeds the precession periods set by the mutual interaction between planets C 
and D, which are of order 10 4 yr. 

If the plane of the stellar binary has a small inclination with respect to the plane of 
the planetary system, we may treat the stellar companion as a third, massive body in our 
linearized equations of motion. It is not difficult to show that the companion star would then 
induce an eccentricity of order (a^/a^eg in the orbit of planet D. For the same parameters 
employed above, this secularly forced eccentricity is of order 10~ 3 . 

Sigurdsson (1992) has shown that the mean induced eccentricity and inclination of a 
planet orbiting a 1.4M star exceed 10~ 3 and 10 -2 , respectively, if a O.7M companion 
approaches the primary to within a distance of ~5 times the semi-major axis of the planet. 
These are mean induced orbital parameters, averaged over 25,000 different possible encounter 
geometries. For the stellar companion having a E ~ 750 AU to come within 5a 2 ~ 13 AU 
of the primary, eg > 0.98. This is an underestimate of the eccentricity required for the 
companion star to be dynamically relevant because m E pa O.2M < O.7M . In any case, 
we consider such large eccentricities to be improbable. Moreover, if such large eccentricities 
were to obtain, the consequent violent forcing of the planetary system over the orbital period 
of the star would threaten the stability of the system. 
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2.2.3. Other Planetary Systems 

We are aware of two other extrasolar two-planet systems whose present-day dynamics 
may be primarily secular: HD168443 (Marcy et al. 2001) and 47 UMa (Fischer et al. 2002). 
In the first system, the less massive, inner planet occupies the more eccentric orbit. This 
is consistent with the idea that both planets excited each other's eccentricities by, e.g., 
divergent resonance crossings (Chiang et al. 2002). Alternatively, each planet's eccentricity 
may have been excited by local disk material (Goldreich & Sari 2002), independently of the 
presence of the other planet. In HD168443, the separation in semi-major axis between the 
two planets is greater than that between planets C and D in v And, increasing the likelihood 
that disk gas was confined between the two planets in the former system. In the absence 
of an innermost disk, ring confinement is necessary for the resonant excitation of the inner 
planet's eccentricity, either by divergent resonance crossings or local disk-planet interaction. 
Both scenarios are consistent with the fact that HD 168443 does not exhibit an apsidal lock; 
there is no reason to expect the ratio e outer / e inner to be driven to large values. 

In the case of 47 UMa, the eccentricity of the less massive, outer planet is constrained by 
observation to be less than 0.2 (Fischer et al. 2002). If said eccentricity is less than 0.1, we find 
that the apsides of the two planets are locked about 0°, consistent with the findings of Fischer 
et al. (2002). Can we apply our theory for v And to 47 UMa? The eccentricities in 47 UMa 
are so low that it is not clear how the system was driven by an external influence, if at all. 
Moreover, the phase space of possible trajectories in Ad)-e outer /ei„ ner space is qualitatively 
different from the case of v And, because most of the angular momentum of the system is 
carried by the inner planet. This causes an island of libration about Auj = 180° to appear 
at e outer I dinner ^ 1- An important feature of 47 UMa is the proximity of this system to a 5:2 
mean-motion resonance; if future Doppler velocity data place the system in this resonance, 
then the apsides circulate (Fischer et al. 2002). Whether or not the system inhabits this 
mean-motion resonance, it seems that our theory for v And cannot be applied to 47 UMa in 
a straightforward manner. 

3. Summary 

We have sketched an evolutionary scenario that explains the origin of the eccentricities 
and the apsidal alignment exhibited by the orbits of planets C and D in the extrasolar 
planetary system, Upsilon Andromedae. In our picture, the eccentricity of the outermost 
planet, D, grows smoothly via resonant interaction with a circumstellar disk of gas remaining 
from the formation of planets C and D. Provided the external, disk-driven amplification 
of the eccentricity of planet D occurs over a timescale that is long compared to apsidal 
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precession timescales set by the planets' mutual gravity, the orbits of planets C and D evolve 
adiabatically through a series of classical Laplace-Lagrange solutions. Increasing the ratio of 
the eccentricity of planet D to that of planet C eventually forces apsidal circulation to give 
way to apsidal libration. Continued pumping of this eccentricity ratio during the phase of the 
libration when this ratio is smallest drives the libration amplitude to zero. The resonantly 
excited eccentricity of planet D is secularly shared with that of planet C. In short, the 
substantial eccentricity of planet C and the apsidal resonance are both consequences of the 
growth of the eccentricity of planet D which, in turn, was driven by a third agent, most likely 
a remnant circumstellar disk lying exterior to the orbits of both planets. The companion 
M4.5V star to v And recently detected by 2MASS is likely dynamically decoupled from the 
planetary system. 

Our mechanism for capture into apsidal resonance and damping of apsidal libration 
is, to our knowledge, novel. Moreover, it enjoys complete independence from energy dis- 
sipation mechanisms that are normally required to damp libration amplitudes as in, e.g., 
the case of the Laplace resonance inhabited by the Galilean satellites (see, e.g., Murray & 
Dermott 1999). The gravitational interactions between planets C and D are, to excellent 
approximation, secular and therefore energy-conserving. 

It is a pleasure to acknowledge Geoff Marcy and Debra Fischer for helpful and friendly 
exchanges and for providing orbital parameters, and Patrick Lowrance and Chas Beichman 
for generously sharing their most recent discovery. 
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